%双扭线
% x=0:0.01*pi:3*pi;
% r=(sin(x./3)).^3;
% x1=3*pi:0.01*pi:4.5*pi;
% r1=(sin(x1./3)).^3;
% x2=1.5*pi:0.01*pi:3*pi;
% x2=x2+3*pi;
% r2=(sin(x2./3)).^3;
% 
% figure(1);
% polarplot(x1,r1);
% hold on;
% polarplot(x2,r2);
% title('r= sin ^3 \theta /3 极坐标')
% 
% figure(2);
% plot(x,r,'b');
% title('r=sin ^3 \theta /3 直角坐标')
% 
% figure(3);
% polarplot(x,r);
x=0:0.01*pi:2*pi;
r=1+cos(x);

figure(1);
plot(x,r);

figure(2);
polarplot(x,r);